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Calculating the Wigner distribution function in the reaction plane, we are able to probe the 
phase-space behavior in time-dependent Hartree-Fock during a heavy-ion collision. We compare 
the Wigner distribution function with the smoothed Husimi distribution function. Observables 
are defined to give a quantitative measure for local and global equilibration. We present different 
reaction scenarios by analyzing central and non- central ^^O+^^O and ^^Zr+^^^Zn collisions. It 
is shown that the initial phase-space volumes of the fragments barely merge. The mean values of 
the observables are conserved in fusion reactions and indicate a "memory effect" in time-dependent 
Hartree-Fock. We observe strong dissipation but no evidence for complete equilibration. 

PACS numbers: 21.60.-n,21.60. Jz 



I. INTRODUCTION 

The time-dependent Hartree-Fock (TDHF) method 
was originally proposed as early as 1930 by Dirac [l[. For 
a long time, it was merely a formal tool to derive many- 
body approaches like, e.g., in Q to derive linear response 
theory. The enormous progress of computational facili- 
ties has made TDHF a practical scheme for describing the 
dynamics of many-body systems. By now it has found 
widespread applications in various areas of physics. Un- 
der the label of time-dependent density-functional the- 
ory it is used in electronic systems like atoms, molecules, 
clusters, and solids, see e.g. The earliest prac- 

tical applications probably appeared in nuclear physics 
[il, where TDHF is a powerful microscopic approach to 
simulate various dynamical scenarios in the regime of 
large-amplitude collective motion, like fusion excitation 
functions, fission, deep-inelastic scattering^nd collective 
excitations; for early reviews see, e.g., [3-[9|. These pio- 
neering applications were still hampered by the computa- 
tional limitations of their time. With the ongoing growth 
of computational power, fully three-dimensional TDHF 
calculations without any symmetry restriction became 
feasible and so renewed the interest in nuclear TDHF, 
for a few recent examples of state-of-the-art TDHF cal- 
culations in many different processes see [iqI-[15|. 

The TDHF approach allows the self-consistent 
quantum- mechanical description of nuclear dynamics on 
a mean-field level. Self-consistency means an unprej- 
udiced description once a reliable energy functional is 
given. This explains the versatility of TDHF. It remains, 
however, an approximation since it is a mean- field the- 
ory. TDHF misses dynamical correlation effects stem- 
ming from nucleon-nucleon collisions, which contribute 
to (two-body) dissipation and thermalization. Their in- 
clusion in a fully quantum mechanical treatment has so 
far only been achieved in homogeneous systems like, e.g., 
[l6l , [l7| . Including dynamical correlations for finite nu- 
clei is presently still restricted to a semiclassical descrip- 
tion |18l-[2Q|. On the other hand, it was found that nu- 



clear TDHF calculations already include a great deal of 
(one-body) dissipation if all terms of the functional, par- 
ticularly the spin-orbit terms, are properly accounted 
for [2l| and if all symmetry restrictions are removed 
(22} . This dissipation within TDHF does not result from 
two-particle collisions but from collision of one particle 
with the boundaries of the moving mean-field potential 
("single-particle dissipation" [23|) which randomizes the 
single-particle states. In a heavy-ion collision, two pic- 
tures of single-particle dissipation can be distinguished. 
The "window" picture describes dissipation of relative 
momentum via nucleon exchange through a neck while 
the "wall" picture deals with the dissipation of kinetic 
energy by reflection of the nucleons at a moving wall 
[25h27J. The latter results in a net increase of the nu- 
cleons' thermal energy provided there is no correlation 
between the nucleonic and wall motions. However, these 
are idealized concepts which are not always immediately 
applicable to realistic heavy- ion collisions (28|-[3Q|. Until 
now it is not understood at a detailed level how rapidly 
and how strongly equilibration works within the TDHF 
approach. 

A rough global measure of dissipation is given by com- 
paring initial and final kinetic energies of the fragments in 
a heavy- ion collision [22|. More detailed analysis should 
look at something like a local momentum distribution. 
This naturally leads to the concept of a Wigner func- 
tion which provides a phase-space picture of a quantum 
state. Originally introduced in j3l|, it is often used for 
establishing the connection between quantum and clas- 
sical physics [32I. The result of such semiclassical limits 
is a mean-field dynamics in classical phase-space called 
the Vlasov equation Q which is widely used in simu- 
lating nuclear dynamics p!8l-[2Q|. In this paper, we want 
to stay at the fully quantum-mechanical level and em- 
ploy the Wigner function as a useful observable helping 
to analyze TDHF dynamics. An early analysis of that 
kind is found in |2^] . The Wigner function has the weak- 
ness that it is not positive semidefinite, thus preventing 
a strict probabilistic interpretation. This defect is cured 
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by some phase- space smoothing leading to the Husimi 
function jsH, [36| , which also turns out to be the better 
starting point for the semiclassical expansion [33|. We 
will also briefly address the Husimi function in connec- 
tion with TDHF results. As the Wigner function is six- 
dimensional and thus rather difficult to handle, we de- 
duce from it more compact measures of dissipation and 
equilibration by considering local quantities integrated 
with some weights over momentum space, e.g., the ec- 
centricity of the momentum space distribution. These 
observables are complemented by others computed with- 
out recurring to the Wigner picture, e.g., the intrinsic ex- 
citation energy which is computed from the local kinetic 
energy density. We will explore these different analyzing 
tools for two realistic applications, collision of ^^O+^^O 
and ^^Zr^^^^Sn. 

The paper is organized as follows: Section [TTl describes 
briefly the numerical handling of TDHF used in this 
work. In Section IHII we present the transformation from 
the TDHF wave function to the Wigner and Husimi rep- 
resentations. Results for the ground states in static calcu- 
lations are presented to compare both pictures. Observ- 
ables are defined in Section [iVl to allow a quantitative 
discussion of equilibration. In Section |Vl we show results 
for dynamical calculations with different nuclei, energies, 
and impact parameters. 

For sake of generality the formal considerations of Sec- 
tion mil are presented in n-dimensional coordinate and 
2n-dimensional phase space. The results in this paper are 
obtained in the reduced two-dimensional reaction plane 
(assumed to be the x-z-plane). For clarity we will label 
the number of coordinate dimensions n of the applied 
distribution function / with f'^^^ 



II. FORMAL AND NUMERICAL FRAMEWORK 

The basis of the TDHF description is a set of occu- 
pied single-particle wave functions ipi{r^t) where / labels 
the states. These wave functions are two-component 
spinors. The Skyrme mean- field Hamiltonian is com- 
puted for_given densities and currents in the standard 
manner [3^]. For all calculations reported here we have 
used the Skyrme parametrization SkI3 (43| . 

The TDHF equations are solved on a three- 
dimensional Cartesian coordinate-space grid. Using the 
fast Fourier transformation (FFT) derivatives can be 
evaluated very efficiently in Fourier-space. The mesh 
spacing is dx = d?/ = dz = 1 fm. 

The stationary ground states of the initial systems are 
computed via the damped-gradient iteration algorithm 
[iol . iZlj . The initial state is obtained by placing the 
ground states of the two fragments in a safe distance and 
giving them a boost towards each other. These states 
are then propagated in time by use of a Taylor-series 
expansion of the time-evolution operator [42] where the 
expansion is taken up to sixth order. The actual time 
step is t = 0.2 fm/c. 



III. WIGNER AND HUSIMI DISTRIBUTIONS 

The Wigner function is a transformation of the density 
matrix to a phase-space function. There are various levels 
of density matrices in a many-body systems and accord- 
ingly various Wigner functions. TDHF can be considered 
as describing the dynamics of the one-body density ma- 
trix p(r, r'), neglecting all correlations between the inter- 
acting nucleons above the mean-field level. This is related 
to the one-body Wigner function which is obtained by a 
partial Fourier transform acting on the relative coordi- 
nate s = r — r', i.e. 

f^\r,k,t) = I ^e-^'^V(r-|,r+|,0 , (1) 

p{ry,t) = J2^]{v,mi{r\t). (2) 

I 

Note that these are, in fact, a spin-averaged density ma- 
trix and correspondingly a spin-averaged Wigner func- 
tion. The dimensionality of the transformation is a 
very compact notation and needs some explanation. Of 
course, our TDHF calculations are always 3D. The full 
Wigner function is then a six-dimensional object, obvi- 
ously a bit bulky. Therefore, we often take cuts and 
look at the Wigner transformation in reduced dimen- 
sions. The notation then means that one coordinate, 
e.g. X, is transformed from the pair in the density 

matrix to the pair (x^kx) in the Wigner function. The 
other two coordinates, y and z in the example, are fixed 
at a certain value yo and zq, usually at the center of the 
nucleus yo = and zq = 0. In other words, f^\x^kx) 
denotes p{x^yo^ zq; x' ^yo^ zo;t) transformed in the x di- 
mension. 

A direct interpretation of the Wigner function as a 
phase-space probability distribution is not possible be- 
cause /w is not positive semidefinite. There can arise 
situations where the quantum oscillations lead to nega- 
tive values. The se p roblems are avoided by the Husimi 
distribution [sH, |36| . The Husimi function /H(r, k, is 
obtained by a convolution of the Wigner function with a 
Gaussian ^(r, k) 

f^\rXt) = J dVd^fe'a(r-r^k-kO 

x/^)(r^k^t), (3) 

a(^Vr-r^k-kO = — e"i^e"i^ , (4) 

ArAk = ^ . (5) 

The Gaussian folding averages /w over the minimal 
phase-space cell of volume {27rh)^ and so successfully 
wipes out the negative values. On the other hand, it 
induces some uncertainty which, however, is physical be- 
cause one cannot localize a particle in phase space better 
than within a volume of (27rfi)^. The Husimi folding has 
one free parameter, the folding width. For best resolu- 
tion in both directions it should be chosen close to the 
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width of the wave functions. As a basis for our choice, 
we use here the nuclear harmonic oscihator model with 
frequency and width parameter given as 

41MeV 



AV3 

This yields the estimate 

2A ~ 2m' 
1 



A = 



mui 



A'/^ fm' 



1 



(6) 
(7) 



4Ar2 2 2A1/3 fm^ * 
The choice is somewhat ambiguous for nuclear reactions 
because one could insert the mass number A for the com- 
pound system or the average A of projectiles, or frag- 
ments respectively. However, these are details which 
do not hamper the analysis; a good order-of-magnitude 
guess suffices for the present analysis. 

In a first round, we investigate the distributions for 
nuclear ground states, to understand the basic pat- 
tern and to have a benchmark from a case certainly 
free of excitation. Figure [T] shows slices through static 
one-dimensional Wigner and Husimi distributions of the 
ground states for three nuclei, a light, a medium heavy, 
and a heavy one. The Wigner distributions show marked 
shell oscillations. The Husimi distributions have effi- 
ciently removed these oscillations and represent a smooth 
curve averaged through the Wigner distributions. The 
amplitude of the shell oscillations decreases with increas- 
ing mass number, but very slowly such that smooth 
Wigner functions (resembling classical phase-space dis- 
tributions) are only reached at an order of magnitude 
A ^ 5000 [H, 113 • The Husimi distributions look smooth 
already for the low mass numbers of really existing nu- 
clei. This, however, is achieved at the price of somewhat 
blurring the details due to the folding procedure. This is 
acceptable for the analysis of the distributions as such, 
i.e. in phase space. It may become misleading when 
reducing the distributions to compact observables by in- 
tegrating over phase space or parts thereof, as will be 
done in Section lIVp . The Husimi folding may add an 
offset to such averaged observables. In such a case, the 
integrations suffice to average out the small-scale oscilla- 
tions in the Wigner functions. Therefore, we will in the 
following concentrate our investigations on the use of the 
Wigner function only. 



IV. MORE COMPACT OBSERVABLES 

The Wigner and Husimi distributions are illustrative 
but difficult to handle, being six-dimensional objects. 
They can be looked at in some selected snapshots and 
by taking cuts through the 6D phase space. Observables 
in lower dimensions down to single numbers are necessary 
complements for the analysis of dynamical processes. In 
this section, we will introduce local observables which are 
distributed in 3D coordinate space. They are reduced to 
single-number observables by further spatial integration. 
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FIG. 1. (color online) Comparison between slices through 
the one- dimensional Wigner f^\z^kz — 0) and Husimi 
f^\z,kz — 0) distribution functions for the static ground 
states of ^^O (top), Zr (middle), and ^^^T/i (bottom). 
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A. Observables from the Wigner distribution 

It is a standard procedure in classical non-equilibrium 
statistical physics to discuss dissipation dynamics in 
terms of the local momentum distribution, i.e. the mo- 
mentum distribution at a given space point [45]. The 
basic features of the local momentum distribution can 
be characterized by its moments. The first moment 



(k(r,^))(n) = 



/d-fe/^)(r,k,t) 



(8) 



plays a special role. It characterizes the center of the 
distribution and it is associated with the average local 
flow. The higher moments are taken as variances, i.e. 
relative to the first moment. For the m-th moment this 
reads 



(k(™)(r,t))(„) 



/d"fc(k-(k(r,t))r/^)(r,k,t) 



(9) 



These moments serve as raw material for further reduced 
observables. Note that they depend on the dimensional- 
ity of the Wigner function used in their definition. This 
is communicated by the index (n) in the moments. 

The radial profile of the momentum distribution may 
be characterized by the ratios of moments, in particular 
the m = 4 to m = 2 ratio 



(10) 



(n) 



Reference value is the thermal equilibrium which corre- 
sponds in the high temperature limit to a Maxwellian 
momentum distribution. These "equilibrium" values are 
given in Table [J for various dimensions. 



TABLE I. Analytic values of the ratio i?*^^-* as defined in 
Eq. (p^ for a Gaussian distribution function, depending on 
the spatial dimension n in which the ratio is evaluated. 



dimension n 


-ftgauss 


1 


3 


2 


2 


3 


5/3 



The ratios are plagued by the fact that cold equilib- 
rium distributions are Fermi functions rather than Gaus- 
sians and, more importantly, are significantly smoothed 
by quantum effects. This hampers an analysis at a de- 
tailed level. A more robust signature of equilibration is 
obtained by the deformation of the momentum distribu- 
tion. The leading term is the quadrupole deformation 
which can be characterized by the eccentricity in the re- 
action plane, which reads 



{kl(v,t))-{kl{v,t)) 
{kl{v,t))^{kl{v,t)) 



(11) 



where the dimensionality index has been skipped for sim- 
plicity. The global eccentricity is obtained by spatial in- 
tegration 



dfc,dfc,/^)(r,k,t), 



(12) 
(13) 



with p the local density. 



B. Intrinsic kinetic energy 

Another interesting observable is the intrinsic excita- 
tion energy. Ideally, it is defined as the difference between 
the actual energy and a "cold" reference energy which is 
obtained from a stationary HF calculation constrained 
to reproduce the density p{v) and current j(r) of the ac- 
tual TDHF state [46|, 113 • The cumbersome density con- 
strained calculations can be avoided when evaluating the 
"cold" reference state in Thomas-Fermi approximation. 
This shortcut was used successfully in Gluster physics 
[38| . The so approximated intrinsic kinetic energy reads 

E\nt(f) = E\^in(t) — £^coll,kin(^) " £^TFw(^) , (14) 

i?kin(t) = ^5]y"d3r|V<p.(r,t)|2, (15) 



coll, kin 



{t) = 



2p{v,t) ' 

£'TFw(i) = y TTFw(r,t) 



rTFw(r,t) = — (37r^)^/Xr,t) 
10m 

I (Vp(r,t))^ 
18m p{r,t) 



^5/3 



(16) 
(17) 

(18) 



It quantifies the non-adiabatic and non-collective com- 
ponent of the kinetic energy, roughly corresponding to 
the intrinsic thermal energy. The first ingredient for the 
calculation is the total kinetic energy, £^kin, of the sys- 
tem. The second term, £^coii,kin(0, subtracts the hydro- 
dynamic kinetic energy contained in the collective flow 
j. The third term, E^Fwit)^ subtracts the instantaneous 
kinetic energy of the zero-temperature ground state at 
the given density p(r,t). The evaluation of this ki- 
netic energy density r(r, t) is done in the Thomas- Fermi- 
Weizsacker approximation 



C. An estimate for the fragment distance 

As a simple observable characterizing the geometry of 
a collisional stage, we introduce the distance d{t) between 
the fragments 



d(t) = Kri(t))-(r2(i))|. 



(19) 




The coordinates of the right and left fragment were 
obtained by sphtting the density of the system symmet- 
ricahy into two half spaces and averaging over each half. 
This is an obvious definition for well separated fragments. 
It becomes somewhat ambiguous in the overlap region, 
but still remains a useful indicator of the overall geome- 
try. 



V. RESULTS 

We present TDHF results for different reaction scenar- 
ios, ^^O+^^O collisions head-on and with finite impact 
parameter, and a ^^Zr-\-^^'^Sn collision. The Skyrme 
parametrization SkI3 [43| is used for the calculations. 
We performed test calculations with other Skyrme forces 
and found very similar results. Thus we report the re- 
sults only from this one force. The central collisions were 
computed on a coordinate space mesh with 48 x 24^ grid 
points and the non-central ones with 36 x 24^ points. 



A. ^^O+'^O Collisions 

1. i^O+i^O Central 

First, we analyze a ^^O+^^O collision with a center- 
of-mass energy of Ec.m. = 100 MeV and zero impact pa- 
rameter 6 = fm. Figure [2] shows the one-dimensional 
Wigner distribution f!^\z^ kz^t) at four different stages 
of the collision. Initially {t = O.Ofm/c), there are two 
cold nuclei far apart from each other. They are shifted 
in /c^-direction depending on their initial boost. At the 
intermediate stage {t = 67.8 fm/c), the phase space vol- 
umes of the two fragments seem to merge but are avoid- 
ing each other, i.e. they maintain a division line. This 
is a consequence of the Pauli principle. After a while 
(t = 155.8 fm/c), the phase space volumes start to sepa- 
rate, keeping some contact still for some time. 

The final stage {t = 273.8 fm/c) shows two separate 
fragments again. But here the structure is quite differ- 
ent as compared to the initial state. Both blobs become 



6 



strongly asymmetric. The kz position of the maximal 
peaks (red spots) are not lower than initially. But the 
asymmetry in the distribution extends very much to- 
wards lower kz and also to values of opposite sign. This 
indicates that the average slowdown in relative cm. mo- 
tion in this case is not due to a global downshift of an oth- 
erwise symmetric distribution, but to the strong asym- 
metry reducing significantly the average kz. A possible 
interpretation is that the wave functions maintain their 
initial momentum structure to a large extent but also 
components from the other fragment are mixed in. 
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FIG. 3. (color online) The global eccentricity s(t) obtained 
from the two-dimensional Wigner function (x, kx, kz^t), 
the distance between the fragments d{t), and the internal ki- 
netic energy Eint{t) for a central ^^O+^^O collision with a 
center-of-mass energy of Ec.m. = 100 MeV. 



In a next step, we analyze the time evolution in terms 
of compact (single number) observables. The time evo- 
lution of the intrinsic kinetic energy Eint(t) and of the 
global eccentricity £(t) are shown, together with the frag- 
ment distance d{t)^ in Figure [3l The time of maximum 
overlap (compound stage) is reached at 75 fm/c where 
d{t) has a minimum. Both kinetic observables show a 
pronounced maximum there. As the reaction continues 
the eccentricity is strongly damped and keeps oscillating 
at a low level. This indicates some thermalization. The 
internal energy reaches its maximum at Eint ^ 86 MeV 
and saturates, again with some persisting oscillations, at 
half of the maximal amount. As the potential energy 
plays a huge role in the compound stage, the values for 
the kinetic energies have to be taken with care here. The 
asymptotic values are more directly interpret able. They 
show significant heating (from Eint{t)) and great deal of 
equilibration already within the short time span of the 
simulation. 

To visualize the oscillations in e{t) and Eint{t) we show 
in Figure [H the two-dimensional momentum distribution 
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FIG. 4. (color online) A kx-kz cut of the two-dimensional 
momentum distribution f^\kx: kz^t) plotted at the center of 
the fragment moving finally with negative mean momentum 
in /c^ -direction. The plots are taken at three different times t 
near the final stage of the calculation. Two contour lines are 
plotted to highlight the levels of /^^ at 2 • 10~^ and 4 • 10~^. 
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/w i^x^ kz^t) in the exit channel for the fragment moving 
to the left (with negative (kz)). The shape is asymmet- 
ric and oscillates back and forth. This indicates that 
the largest collective effect in the exit channel is residual 
octupole oscillations which have their counterpart also 
in similar octupole oscillations of the fragments' spatial 
shape. 

We have checked the momentum ratio R^^^ (r, t) as 
given in Eq. (fTO]) for n = 1,2 at different times to 
probe the closeness of the momentum distribution to a 
Maxwellian distribution. The comparison with the ana- 
lytic values i^gauss from Tab. [J is not conclusive, as quan- 
tum effects blur the classical concept behind this ratio 
by making the distributions too different from Gaussians. 
This casts some doubts on the usefulness of the global ra- 
tio R{t) in this still predominantly quantum-mechanical 
domain. We will come back to this observable in Section 

lYll 



2. ^^O+^^O Non-Central 

In this Section we analyze ^^O+^^O fusion reactions 
with non-zero impact parameter 6 = 2 fm and two differ- 
ent center-of-mass energies Ec.m. = 20 MeV and Ec.m. = 
50 MeV. 

Figure [5] shows the global eccentricity e{t)^ the inter- 
nal kinetic energy Eint{t)^ and the fragment distance d{t) 
for both reactions. The distance oscillates after the min- 
imum reached on first impact. This indicates that both 
cases describe a fused compound state. The eccentric- 
ity £(t) follows the oscillations of the distance, reflecting 
a continuing vivid interaction with the spatial deforma- 
tion. This indicates that we are far from equilibration. 
The intrinsic kinetic energy grows initially and soon lev- 
els off, leaving small residual oscillations about a constant 
mean value. This mean intrinsic energy is, of course, 
larger for the higher-energy collision {Ec.m. = 50 MeV). 
The example demonstrates nicely that one needs a cou- 
ple of observables to conclude on equilibration. One may 
be tempted to take the constant Eint{t) as indicator of 
a thermalized state. The still large values of eccentric- 
ity and the oscillations thereof prove clearly that we are 
rather in a situation of substantial coherent oscillations 
of the compound system. In this context, it is to be re- 
membered that the energy stored in the collective motion 
of the compound system is subtracted in the evaluation 

of Eint{t). 

The demonstrated behavior of the observables was 
checked up to t = 4000 fm/c, twice the time span shown 
in Figure [5l The pattern carried on unchanged also for 
these longer times. 

B. 9^Zr+i32g^ 

As an example for a much heavier nuclear system we 
present fusion of ^^Zr+^^^Sn achieved with a center-of- 
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FIG. 5. (color online) The global eccentricity e{t) obtained 
from the two-dimensional Wigner function (x, z,kx,kz,t), 
the distance between the fragments d{t), and the internal ki- 
netic energy Eint{t) for two "^^0+"^^0 fusion reaction with a 
center-of-mass energy of Ec.m. = 20 MeV (top) and Ec.m. = 
50 MeV (bottom) and impact parameter h — 2 fm. 



mass energy of Ec.m. = 250 MeV and impact parameter 
h = 2 fm. The reaction between the neutron-rich ^^^Sn 
nucleus and ^^Zr was already studied in TDHF with a 
focus on barrier heights and widths of the heavy-ion po- 
tential as well as capture cross sections (4^ . 

It was not possible in the present analysis to calculate 
the distance d{t) between the fragments as it was done in 
the ^^O+^^O fusion scenario. The numerical algorithm 
selecting the spatial expectation values for a two-body 
system was not able to detect two distinct objects dur- 
ing the whole calculation and in this asymmetric system 
a simple symmetric division of the grid was not possi- 
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FIG. 6. (color online) The global eccentricity €{t) obtained 
from the two-dimensional Wigner function f^-* {x, z^kx^kz^t), 
the quadrupole Q20, and the internal kinetic energy Eint{t) 
for a ^^Zr+^^^Sn fusion reaction with a center-of-mass energy 



of Ec 



250 MeV and impact parameter 6 = 2 fm. 



ble. In Figure [6] we therefore use the expectation value 
Q20 = {Q20) of the quadrupole operator Q20 to visualize 
the global geometry of the reaction. Again, large values 
indicate separated fragments and low values a compound 
stage. It is obvious from the figure that the reaction ends 
in a compound nucleus. The overall trends of intrinsic 
kinetic energy and eccentricity are to some extent sim- 
ilar to the results in Figure [5l However, the final ec- 
centricity is much smaller, still maintaining some small 
oscillations. This indicates a better thermalization than 
seen for ^^O+^^O, which is no surprise because the single- 
particle phase space is much larger for the heavier system. 
The trend of the intrinsic energy does also differ in de- 
tail. There seem to be two stages of growth, a fast initial 
rise on the way to the compound stage and a slower, but 
steady, growth up to 1000 fm/c. This indicates that some 
thermalization processes and energy transport from de- 
formation energy to kinetic energy is still going on. After 
1000 fm/c we again see a rather constant Ei^t as seems 
to be typical for energetic compound nuclei. 

Figure [71 shows the local ratio R^'^\r^t) as defined 
in Eq. (p!Q|) in the reaction plane at t = 1279.8 fm/c. 
The surface region is distinguished by large values com- 
ing close to the Maxwellian reference values while much 
smaller ratios are seen inside. In order to illuminate these 
results, we have a closer look at the more detailed mo- 
mentum distributions at four selected points indicated in 
Figure H Figure [8] shows results for two the outer points 
at the surface. The first point (Pi) is taken at approxi- 
mately half the maximum value of /w- The distribution 
is very similar to a Gaussian overlayed by a slight asym- 
metry. The next point (P2) shows a more pronounced 



FIG. 7. (color online) The local ratio R'''^\r) obtained from 
the two-dimensional Wigner distribution f^\x, z^kx^kz^t) 
for a ^^Zr+^^'^Sn fusion reaction at t = 1279.8 fm/c. Four 
points are marked in this plot to be analyzed later more pre- 
cisely (Figure [8] and Figure [9|). The reference value from a 
Gaussian distribution is i?giuss = 2. A Contour line is plotted 
to highlight the level of R^^\r) = 1.5. 



asymmetric shape. Moving further to the inner points 
(P3,P4) reviewed in Figure [9] the momentum distribu- 
tions differ substantially from Gaussians and come closer 
to the idea of a Fermi distribution, although heavily over- 
layed by quantum shell oscillations. 

An similar analysis of the momentum distribution at 
other points near the nuclear center yields similar results. 
The inner region of the merged systems seems to stay 
rather "cold" during the reaction. 

The strong quantum mechanical shell oscillations hin- 
der a fit of the distribution functions shown in Figure 
[8] and [9] to a Fermi function from which one eventually 
could read off an estimate for the system's temperature 
distribution. Therefore we compare the ^^Zr+^^^Sn sys- 
tem to be assumed "hot" with the "cold" analogue of 
this system. Figure [10] shows the momentum distribu- 
tion of the ^^Zr-h''^^^Sn compound system at the point P4 
indicated in Figure [9l This is compared with the result 
from the prolate ground state of ^^OTh. The ^^Zr+^^^g^ 
system consists of p = 90 protons and n = 138 neutrons. 
230 rpj^ shares the same proton number with two addi- 
tional neutrons. The ground state nucleus shows huge, 
fully developed shell oscillations. Compared to these, 
the remaining quantum oscillations in the "hot" com- 
pound state become rather small. The disappearance 
of quantum shell effects is a major thermalization effect 
[32l,l48|. The occupation of high momentum components, 
however, which would also be expected for hot systems 
remains insignificant. This is due to the fact that the nu- 
cleus is an open system from which high energy particles 
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FIG. 8. (color online) The left column reviews the two-dimensional momentum distribution f^\kx, kz) in the selected (outer) 
points from Figure [71 Contour lines are plotted to highlight the levels of /^^ at 1 • lO"'^, 2 • 10~^, and 3 • 10"^. Slices through 
f^\kx, kz) matching the kx, kz-a^xis are shown in the right column. 



escape, constantly depleting the high-momentum parts 
of the distribution. This explains why a fit to Fermi dis- 
tributions failed. A measure of temperature may be de- 
duced from the suppression of the shell oscillations, but 
this analysis is blurred by the large thermal fluctuations 
in the momentum-space density. For the time being, the 
eccentricity remains the cleanest indicator of equilibra- 
tion. 



VI. SUMMARY 

In this work we have analyzed from different perspec- 
tives the dynamics of TDHF during various reactions in- 
cluding the nuclei ^^O, ^^Zr, ^^'^Sn, and '^^^Th with vari- 



ous center-of-mass energies and impact parameters. The 
key quantity of the analysis is the Wigner distribution 
function which provides a detailed phase-space picture 
of the quantum state. As complementing quantities, we 
also considered three more compact observables in terms 
of local distributions: the ratio i?(r, t) of the weighted 
moments (weight four and two) of the local momentum 
distribution described by the Wigner function (i.e. in- 
tegrating the Wigner function over momentum space for 
fixed local position), the eccentricity e(r, t) of the local 
momentum distribution, and the intrinsic excitation en- 
ergy £^intr(r, t) as deduced from the kinetic energy den- 
sity. 

General properties of the Wigner distribution were dis- 
cussed first for stationary states. It shows oscillations 
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FIG. 9. (color online) Same as Figure [8] but for the selected (inner) points in Figure [71 
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which stem from the quantum shell oscillations of the 
underlying single-particle states. We also looked at the 
Husimi function derived from the Wigner function by 
some phase-space smoothing. The latter indeed provides 
a cleaner and more intuitive picture. We find, however, 
that the shell oscillations are much reduced in the dy- 
namical scenarios of heavy- ion collisions, allowing us to 
continue the dynamical studies with the Wigner function 
alone. 

We have visualized the collision process through snap- 
shots of a 2D cut of the 6D Wigner function. This shows 
that the two initially separated phase space blobs never 
fully merge, even at the compound stage. The distribu- 
tions of the emerging fragments acquire a strong asym- 
metry in momentum space and nicely show the phase 
space rotations associated with the remaining octupole 
oscillations of the final fragments. 



The moment ratio R was intended as a means to com- 
pare the shape of the TDHF- Wigner distribution with a 
Maxwellian distribution corresponding to thermal equi- 
librium. We find for all reaction parameters that the 
moment ratio remains below the Maxwellian reference 
value, which means that thermalization could not be as- 
serted in this observable. The reason is that that the 
high-momentum tail of the actual distribution is imme- 
diately depleted by particle emission. This exemplifies 
the fact that a true equilibrium state is hard to establish 
in an open system. 

The eccentricity e turned out to be a more useful in- 
dicator. It grows dramatically in the initial phase of the 
reaction and relaxes to lower values quickly after the com- 
pound state and then remains oscillating about some fi- 
nite value. This means that the final relaxation to a 
thermal state is probably underestimated in mere TDHF. 
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FIG. 10. (color online) Slices along kx (left) and kz (right) through fy^ (kx, kz) are shown for the merged Zr-\-^^^ Sn system 
(red) at t = 1279.8 fm/c in the point P4 from Figure [Tlin comparison with the ground state of ^^^Th (blue). 



Similar patterns are shown by the intrinsic kinetic energy 
Eint' The resulting "asymptotic" value of Eint depends 
strongly on the initial conditions, e.g., growing with the 
initial collision energy. 

We conclude that although TDHF includes dissipation 
owing to single-particle viscosity, which acts strongly in 
the initial phase of reactions, there is no evidence for 
complete equilibration. 
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